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We discuss the standard approach to the problem of the low momentum limit of the spectrum for 
a weakly interacting Bose gas. The Bogoliubovs spectrum is shown to be obtained as a Goldstone 
mode thanks to the introduction of a chemical potential \i. This procedure has, however, difficulties 
since the breaking of the gauge symmetry implies that the corresponding chemical potential must 
be taken as zero, unless it is introduced before breaking the symmetry. But if this is done, after the 
symmetry breaking /1 loses its meaning as a chemical potential. An alternative two-mode solution is 
• suggested having two modes, one of them being the free-particle quadratic in momentum spectrum, 

the second bearing a gap. This gap leads to a A-type behavior of the specific heat near the critical 
temperature. 
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I. INTRODUCTION 
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Most studies for the weakly interacting Bose gas deal with the zero temperature approximation. We must cite at 
first the famous Bogoliubov model . In it the assumption N — n = ^2 p n p <C N is made, so that one can discard a 
term of order {J2 P n p) 2 - This leads to a gapless spectrum, but such approximation is actually valid for temperatures 
extremely close to 0°K. Bogoliubov's sprectrum is in agreement with the pure phenomenological predictions made by 
Landau concerning the superfluid spectrum, Q whose lower branch is expected to have a phonon behavior. Far from 
0°K, as happens near and below the critical temperature T c , the approximation made by Bogoliubov is, however, not 
qq , valid, and the spectrum bears a gap |3(, since it is, 

o ■ , 

\Q ; e(p) = y/{p 2 /2m) 2 + 8K(p 2 /2m) + 12K 2 , (1) 

o 

where K = Uono/2V, and Uo is the (approximately) constant potential, no the condensate and V the volume. As 
shown in this gap, not being satisfactory with the expected spectrum, bears however the interesting property of 
leading to a divergent behavior of the heat capacity of the system, which fact is in correspondence with the well-known 
lambda- type phase transition of He 4 . The usual Bogoliubov gapless spectrum does not predict any A behavior |4| and 
is valid only near T = 0, also it does not exhibit even the discontinuity in the C v derivative at T c typical of the usual 
q I ideal Bose gas. 

Bogoliubov's work was a brilliant insight in the appearance of gapless particles after spontaneous breaking of a 
gauge symmetry, a fact which became more clear 15 years later, after the Goldstone theorem |15|. After Bogoliubov 
the weakly interacting gas spectrum was investigated by Beliaev by using the Green's function method, but the 
^ correspondence with Bogoliubov's results could be obtained only after introducing a quantity identified as a chemical 
potential, which was not taken into account in the original model by Bogoliubov. Below we shall see, by a different 
method than the one used by Beliaev, that spontaneous symmetry breaking is produced by introducing explicitly 
the chemical potential, leading to the appearance of a gapless Goldstone mode. On the other hand, by taking the 
chemical potential as equal to zero in Beliaev's calculations, the spectrum Q is reproduced. 

In the present letter we are interested essentially in discussing the meaning of taking such a chemical potential as 
different from zero, at the light of the fundamental principles of statistical mechanics. 
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II. THE HAMILTONIAN 



Condensation breaks the gauge symmetry (sec below). The introduction of a chemical potential coupled to a 
Noether charge Q associated to a broken symmetry in the density matrix was investigated in the relativistic context 
in papers -8|. Such charge has many odd features from which we mention two: first, in general, it is not conserved, and 
second, its introduction ban the fulfillment of the Kubo-Martin-Schwinger (KMS) condition required for the statistical 
description of systems in thermodynamical equilibrium, since according to the Araki-Haag-Kastler-Takesaki (AHKT) 
theorem Q , only those charges that belong to the center of the unbroken part of the symmetry group satisfy the KMS 
condition and can have a chemical potential associated to them. 

In the non-relativistic weakly interacting Bose gas case, the charge operator associated to the broken symmetry 
is N — J d 3 xip*(x)4>(x) — J2p=o° a *(p) a (p)i where it is taken the condensate density as a*(0)a(0) = no and in 
consequence a(0) = v / n^e ie ,a*(0) = y / n^e _ie .This means ^(J that the physically relevant representations of the 
algebra of canonical variables are labelled by the parameters no, 0, and 

limy^oo J d 3 x < no, 9\tp(x)\n ,8 >= y/n~o~e ie (2) 

which for n ^ is not invariant under gauge transformations 
Let us consider the Lagrangian 

C = -r(x)id Tp(x) + r(x)^(x) - \(r(xMx)) 2 (3) 

2m 

where A = |^ and we are considering the approximation of a constant repulsive potential. Under the gauge transfor- 
mation 



if>(x) -» e ia ip(x) (4) 

terms having equal number of ip(x) and ip*(x) are gauge invariant, so C is gauge invariant, as it is also the density of 
particles N. 

From the Lagrangian we can get immediately the Hamiltonian in which, by taking the Schrodinger representation 
of the fields ip(x) = y ~}2 a(p)e^' r ^ h , ip*{x) — y ^ a*(p)e~ l P' r / h , one obtains in momentum space an expression 
quadratic in the fields a*(p) , a(p) . After condensation, if it is taken a*(p) = a*(p)e~ t6 and a(p) = a(p)e obviously, 
such a Hamiltonian is not gauge invariant. If we fix the gauge 9 = 0, we get finally the bilinear part of the Bogoliubov's 
Hamiltonian, 
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H B = Xnl + y\a*(p)^-a{p)+K{Aa*{p)a{p)+a*{p)a*{-p)+a{p)a{-p)] (5) 
t-^ 2m 

v 

where K = Xtiq. Note that the coefficient of the first "potential" term is 4K and not 2K, as is done as an approximation 
in the usual treatment of the problem in terms of the total number of particles N. 



III. DENSITY MATRIX AND BELIAEV EQUATIONS 



The immediate task is to build a density matrix aimed to obtain the partition functional Z. If we name N p = 
J2 P a *(p) a (p) the operator for the number of excited particles, we see that [H B ,N p ] ^ 0. In consequence, there is no 
a common set of eigenstates for these two operators and the density matrix p = er^ 13 ^^^ 13 is an ill-defined one for 
describing states of equilibrium. If we define fl = —ft^ 1 InZ one cannot state that N = — -^f, is the average number 
of excited particles. However, one might argue that we have the right to consider /i ^ as a Lagrange multiplier, 
since in building p and O it loss the meaning of a chemicalpotential but then we face new difficulties. 

Working according to Bogoliubov's standard procedure 0, we may proceed to diagonalize Hb, 

H B =J2<P) b *(P) b (p) (6) 
p 

where according to Bogoliubov's transformation, we have b*(p) = [a(p) + a p a*(—p)]/(l — a 2 (p)), b(p) = [a*(p) + 
a p a{—p)]/{\ — a 2 (p)), and a(p) = [AK + p 2 /2m — e(p)]/2K). In this representation, in which b*(p)b(p) is gauge 
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invariant under the transformation b(p) — > e l9 b(p), b*(p) — > e l9 b*(p), we have 

(1 + a(p) 2 )n p + a(p) 2 a % 



N = J2- 



1 — a(p) 2 1 — a(p) 



(6*(p)6*(-p) + 6(p)6(-p)) J (7) 



which is obviously not gauge invariant. The conditions of the AHKT theorem Q are not satisfied. Thus, the system 
of equations l|17|l and l|18|l are also not gauge invariant. 

Thus if we continue in a pure formal way, and leave /i ^ 0, by considering Hb — fJ>(N p + n ) as a new Hamiltonian, 
one can obtain some partition function Z from it, from which the system of Green functions and the thermodynamical 
potential can be obtained. By introducing adequate external sources, we have for the partition function obtained 
from the unperturbed Hamiltonian Hb, 



2(7], C) = M(P)e- f,VXr *-' ma Zp (8) 
where M(/3) is a temperature dependent constant and 

Zp = J Va*Vae^J S(p)+n{p)a ' ip)+C{p)a ^ (9) 

and 

S{p) = (ip 4 + ^-^)a*(p)a(p)-K(4a*{p)a{p) + a*{p)a*(-p)+a(p)a(-p)) (10) 

We may consider the interaction term in the Hamiltonian, containing terms of third and fourth order in a*(p),a(p). 

If V(a(0),a*(0),a*(p),a(p)) is such term, the interacting partition function is Zj = e ^Pi.p V{a{0),a ' {0),J ^ J ^ ) Z. By 
using the standard methods of functional differentiation, we may develop a formalism to obtain from Zj the equations 
for the Green functions G(p,p') — 5 2 Zi /8((p)5r](p r ), G\(p,p') = 8 2 Zi / 8rj(— p)8r)(p'). We will restrict, however, to the 
corresponding expressions for the V — approximation Zj — Z, by taking the bilinear terms in l |10| l. and by calling 

Gq 1 = (ip 4 + M — 2m) one easily gets, 

G(p) = G (p) + 4KG (p)G(p) + 2KG a (p)G 1 (p) (11) 



dip) = 2KG (~p)G(p) + 4JTGf (-p)Gfi(p) (12) 
These are a simplified version of the Beliaev's equations 5|, which lead to 

g 1 ( P ) = ?fy d3) 

v ' G^i-p)-^ v ; 

( Pi -ie(p)){p 4 + ie(p)) 

where 



e( P ) = \l(^-fL + 4K) 2 -4K 2 . (15) 
V 2m 

The effective potential is obtained as U = /3 _1 InZ/V, by defining ^ = y/no, as 

U = A£ 4 - M C 2 + 0/3 (16) 

where 0/3= \nZp. We see that the chemical potential /i plays an essential role in showing explicitly the symmetry 
breaking term in U, analogous to the "negative mass term" of the usual relativistic theory [15|]. Its presence leads 
necessarily to the existence of Goldstone bosons. 
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From l|16(l one can obtain the condition of extremum of U 



d -g = 4A{= - 2K + ^ = 0, (17) 



But we cannot write that the quantity 



d^ = no + ^T (18) 



as the average number of particles, i.e, N q = — ^Sf- is not the number of excited particles. Usually it is interpreted as 
the average number of " quasiparticles" . For T — ► 0,N q — ► and only in that case no ~ N. 

The equation (|17f) gives £ 2 = /Lt/2A as the minimum at the tree level, i.e., for T — [i = 2K. By substituting in l|15|l . 



this leads immediately to the (gapless) Goldstone boson c(j>)b — y (§^) 2 + Jm' which is the usual Bogoliubov 
spectrum. According to Goldstone theorem, these modes if calculated correctly, would maintain gapless at all orders 
of perturbation theory. This result is in correspondence to the one obtained, before the formulation of the Goldstone 
theorem, by Hugenholtz and Pines [Tl| . 

We see that for the reasons pointed out above, there is no right to build a density matrix in the grand canonical 
ensemble by subtracting from Hb the operator N p multiplied by a nonzero chemical potential /i, since as [Hb, N p ] =/= 0, 
there is no a common set of eigenstates and eigenvalues for them, and p = e - ( ^fB- ' lA p' ,3 does not satisfy the KMS 
condition. Thus, it is unsatisfactory to attribute a standard physical meaning to quantities derived from it if fi ^ 0, 
in particular, to interpret in i|18|) as the average number of excited particles. We introduce the notion of 
"quasiparticles" in part to avoid such a difficulty. We point out also that if we take fi = (|15|) turns into (1). 



IV. ALTERNATIVE PROCEDURE 



We conclude that the standard introduction of a chemical potential is not satisfactory if it is done after breaking 
the gauge symmetry. If we include it before the symmetry breaking, it loss its meaning as a chemical potential 
and becomes a parameter essential for the fulfillment of the Goldstone theorem. However it leads (according to the 
gau ge fixing) to other infinite non-equivalent set of representations of the operator algebra, and of their dynamics 
|l0| . There is no a priori criteria to state which of these representations is physically better. One may choose, for 
instance, that one wants the energy eigenvalue E(p) to be linear in p for p — > 0. Then one possible solution is the 
Bogoliubov-Beliaev spectrum. But from the point of view of Goldstone theorem and its consequences, the use of the 
fields ip*( x ) an d ^{ x ), both containing the symmetry-breaking parameter, is rather obscure. Only one physical mode 
is present. It would ban the appearance of a pure gapless mode and a mode bearing a gap, which is to be expected 
(fron the Goldstone theorem point of view) since we are using two-component fields. 

If we do not tie to the linear mode requirement, an alternative model is, for instance, the one obtained by taking 
in |J3} ip = £ + <j(x) + i h(x). At the tree level these two modes have the spectra Eh = p 2 /m and E a = p 2 /m + A a 
where A CT = K = Xno and no is the condensate (which is the total number of particles at T — 0) The first mode has 
the property of being gapless, and trivially corresponds to the Goldstone Boson 0, since Eh — > for p — > 0. The 

second mode bears a gap. At T ^ 0, this gap is expected to be of form A CT (T) = A n ^/l — (^) 3 , where T c is the 

critical temperature for condensation. It bears the interesting property of leading to a A - behavior near T c ,as shown 
below. We observe that the Bogoliubov-Beliaev spectrum is just the geometric average of the spectra E(p) = y/E a Eh 
this spectrum does not lead to an A-type behavior. 

The divergence of the specific heat cy is seen Q easily starting from E a . One can write the thermodynamic 
potential as — AT p 2 dp \n(l — e -6 ™' 3 ), where A is a constant. By taking x = p/\/2mT, in the infrared limit, 
where po is a sufficiently small momentum, we have 

A CP" p 4 dp A{2mTf/ 2 yw x 4 dx A{2mfl 2 T^y 

rj2 Q 1/2 

From c v — — T|^, we get that cy diverges as A CT for T — > T c . 

It may be argued however, that the existence of a gap, as well as the non-linear gapless mode, contradicts the 
generally accepted idea of a phenomenological linear in p phonon spectrum . One concludes that the model we are 
considering is an oversimplification of the real system, an that one must study the higher order corrections. Work 
in this sense is in progress. Another modification would be to introduce some term accounting for the oscillatory 
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motion of the condensed or quasi-condensed atoms (the spectra obtained above refers actually to translational modes), 
leading to a phonon field having linear behavior of the spectrum. New models for superfluidity are currently being 
elaborated [Tj| on much more phenomenological basis by including such vibrational terms. 
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